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Abstract 

We use molecular dynamics computer simulations to study the relaxation 
dynamics of a viscous melt of silica. The coherent and incoherent intermedi- 
ate scattering functions, Ffi{q,t) and Fs{q,t), show a crossover from a nearly 
exponential decay at high temperatures to a two-step relaxation at low tem- 
peratures. Close to the critical temperature of mode-coupling theory (MCT) 
the correlators obey in the a-regime the time temperature superposition prin- 
ciple (TTSP) and show a weak stretching. We determine the wave-vector 
dependence of the stretching parameter and find that for Fd{q, t) it shows os- 
cillations which are in phase with the static structure factor. The temperature 
dependence of the a-relaxation times r shows a crossover from an Arrhenius 
law at low temperatures to a weaker T-dependence at intermediate and high 
temperatures. At the latter temperatures the T-dependence is described well 
by the power law proposed by MCT with the same critical temperature that 
has previously been found for the diffusion constant D and the viscosity. We 
find that the exponent 7 of the power law for r are significantly larger than the 
one for D. The wave-vector dependence of the a-relaxation times for F(i{q,t) 
oscillates around r(g) for Fs{q,t) and is in phase with the structure factor. 
Due to the strong vibrational component of the dynamics at short times the 
TTSP is not valid in the /3-relaxation regime. We show, however, that in this 
time window the shape of the curves is independent of the correlator and is 
given by a functional form proposed by MCT. We find that the value of the 
von Schweidler exponent and the value of 7 for finite q are compatible with 
the expression proposed by MCT. Finally we discuss the g-dependence of the 
critical amplitude and the correction term and find that they are qualitatively 
similar to the ones for simple liquids and the prediction of MCT. We conclude 
that, in the temperature regime where the relaxation times are mesoscopic, 
many aspects of the dynamics of this strong glass former can be rationalized 
very well by MCT. 
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I. INTRODUCTION 



If a glass forming liquid is cooled from high to low temperatures one finds that its 
relaxation time r, or its viscosity rj, increases rapidly by many decades If r]{T) is plotted 
in an Arrhenius plot, i.e. \og{rj) versus 1/T, one finds that the shape of the curves is not 
universal but depends on the material. For two extreme cases of the shape, Angell has coined 
the terms "strong" and "fragile" glass formers P). Strong glass formers are liquids whose 
Ti{T) curve shows an Arrhenius law in the whole accessible temperature range. In contrast 
to this, fragile liquids show a pronounced crossover at intermediate temperatures from an 
Arrhenius-like law at high temperatures to an Arrhenius-like law at low temperatures with 
a higher activation energy. It is one of the merits of the so-called mode-coupling theory 
(MCT) of the glass transition |]^,^ to offer an explanation for the existence of this crossover 
in terms of non-linear dynamical feedback effects and to make detailed predictions for the 
dynamics of the glass forming liquids close to this crossover temperature Tc, the critical 
temperature of MCT. In the past years there have been a large number of attempts to check 
the validity of these predictions and the results of these efforts was that the theory is indeed 
able to give a surprisingly good description of the dynamics in this temperature range [Q. 
It has to be emphasized that some of these tests concerned not only the predictions of the 
theory on a qualitative level, but also on a quantitative one, and that also the outcome of 
these tests were often in very good agreement with the theoretical prediction. Thus one can 
say that the theory is able to describe in the vicinity of the dynamics of fragile liquids on 
a qualitative as well clS cl quantitative level. 

Things are much less clear for the case of strong glass formers since by definition they 
do not show a crossover between two temperature dependences in r]{T), or only a very mild 
one. Hence it is a priori not clear at all whether these types of systems have a crossover 
temperature and hence whether there is a temperature regime in which MCT can be used 
to describe their dynamics. It was therefore a bit surprising when it was found that glycerol, 
a glass former which is a rather strong glass forming liquid, shows a dynamics which can 
be described well by means of this theory |0J^. Very recently a computer simulation study 
of liquid Si02, the paradigm of a strong glass former, showed that the dynamics of this 
system at high temperatures shows features which are very reminiscent of the dynamics of 
fragile systems around their Tc 0. In particular it was shown that also in silica a critical 
temperature Tc can be identified, with a value around Tc = 3330 K, i.e. a temperature 
which is presently outside the reach of experimental investigations. It is of interest to note, 
however, that for this material an extrapolation of experimental data to higher temperatures 
gives evidence that around 3220 K a bend in the viscosity would indeed be observed [§, thus 
supporting the results from the simulation. Very recently it has also been demonstrated that 
MCT can be used to understand the dynamics of liquid silica not only on a qualitative level, 
but also on a quantitative one. E.g. it was shown that the theory is able to predict with good 
accuracy quantities like the wave- vector dependence of the Debye- Waller factor 0. The only 
input in that calculation were static quantities, i.e. the partial structure factors and static 
three-point correlation functions, which were determined directly from simulations. Also 
for water, another network forming liquid, it was found that MCT predicts the wave-vector 
dependence of the Debye- Waller factor with a very good accuracy [|10|. Thus these studies 
give evidence that MCT is also able to give a correct description of the dynamics of strong 
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liquids in the temperature region where the dynamics starts to slow down significanly, hence 
expanding the range of applicability of the theory considerably. 

From the point of view of the theory this possibility is perhaps not that surprising. When 
the theory was developed originally, some terms in the equations of motions were dropped 
in order to simplify the equations and because it was argued that they are not very relevant 



to understand the main predictions of the theory This "ideal version" of the theory has 



subsequently been tested extensively in experiments and computer simulations A bit 

later it was recognized, however, that close to the critical temperature the terms that were 
dropped will change some of the predictions of the theory considerably [^,|T^], since they 
are needed to describe processes, today called "hopping processes", that become relevant 
at low temperatures. If these hopping processes are very small, the predicted temperature 
dependence of the relaxation times shows a strong bend in an Arrhenius plot, reminiscent 



to the one seen in fragile glass formers ||13[. If these processes are very pronounced this 
temperature dependence shows only a weak bending around T^, thus a similar behavior as 
the one found in strong glass formers. Thus it seems that a priori the theory is indeed able 
to describe also strong glass formers. 

However, it is presently not clear whether the extended form of the theory, i.e. the ver- 
sion of the theory which tries to include hopping effects, is really reliable also beyond a 
phenomenological description. Therefore it is important to check to what extent the the- 
oretical predictions survive in the case where one deals with a real material or a realistic 
microscopic model in which hopping effects are expected to be important. Such investi- 
gations might also give an idea how the theory can be improved or extended in order to 
describe the dynamics of realistic glassforming liquids. The present paper is hence a contri- 
bution to this issue. For this we use molecular dynamics computer simulations to investigate 
the dynamics of liquid silica in great detail. In particular we will study to what extent the 
behavior predicted by the theory can be found in this dynamics, i.e. we test in a stringent 
way to what extent MCT is able to give a reliable description of the relaxation dynamics of 
strong glass formers at high temperatures. 

The rest of the paper is organized as follows. In the next section we give the details on 
the model we use and the simulations. The following section is devoted to the presentation 
of the results, which are summarized and discussed in the last section. 



II. MODEL AND DETAILS OF THE SIMULATION 

In this section we present the model we used to describe the interaction of silica and 
give some of the details of the simulation. Further information on the latter can be found 



in Ref. 14 



In order to obtain a realistic description of the structural and dynamical properties of 
a given material by means of classical simulations, it is necessary to have a potential at 
hand which describes the interactions between the different ions realistically. Due to the 
importance of silica in applied as well as fundamental science it is nowadays possible to find 
a multitude of different potentials in the literature, most of them having certain advantages 
and disadvantages. Various simulations have shown that one of the most reliable of these 
potentials is the one proposed by van Beest, Kramer, and van Santen (BKS) which was 
obtained by determining the energy surface of a Si04H4 tetrahedron by means of ab initio 
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calculations, parametrizing this surface with a simple analytical expression and subsequently 
doing some lattice dynamics calculations with it to improve the fit parameters [|T5[] . 
The functional form of the BKS potential is 

Mr) = ^^ + A^pexp{-B^^r)-^ a,f3e[Si,0], (1) 

where r is the distance between the ions of type a and (3. The values of the constants 
Qa, qi3, Aa/3, Ba/3, and Ca/3 cau be found in Ref. [l^. For the sake of computational efficiency 



the short range part of the potential was truncated and shifted at 5.5 A. This truncation 
has also the benefit to improve the agreement between the density of the amorphous glass at 
low temperatures as determined from the simulation with the experimental value. Previous 
simulations have shown that this potential is able to reproduce a variety of properties of 
real amorphous silica, such as the density anomaly, the static structure factor, the specific 
heat, the viscosity and the thermal conductivity [0. Thus it can be expected that it is also 
sufficiently accurate to give a reliable description of the time dependence of the intermediate 
scattering function, the central quantity discussed in the present paper. 

The system investigated has 8016 ions in a cubic box with fixed size L = 48.37 A. 
This size is sufficiently large to avoid the finite size effects found in the dynamics of such 
systems [p!7| , p!8| . The equations of motion have been integrated by means of the velocity form 
of the Verlet algorithm, using a time step of 1.6 fs. The Coulombic part of the potential 
has been evaluated by means of the Ewald sum using a constant aL=12.82. Before the 
microcanonical production runs were started the system was equilibrated by coupling it to 
a stochastic heat bath. The duration of this equilibration was significantly longer than the 
typical relaxation times of the system. Therefore we are sure that all the results presented 
below refiect the relaxation dynamics in equilibrium, i.e. that they are not affected by aging 
phenomena. The temperatures investigated were 6100 K, 5200 K, 4700 K, 4300 K, 4000 K, 
3760 K, 3580 K, 3400 K, 3250 K, 3100 K, 3000 K, 2900 K, and 2750 K. In order to improve 
the statistics of the results we averaged at each temperature over two independent runs. 

Note that all these temperatures are above the melting temperature of real silica, which is 



around 2000 K [|I9| , a temperature which the BKS potential is able to reproduce reasonably 
well [^]. Thus none of our simulations probe the dynamics of the system in its supercooled 
state. Nevertheless we will see below that even at these relatively high temperatures the 
dynamics of the system is very slow. Thus this shows that in order to have a slow relaxation 
dynamics it is not necessary to be in a supercooled state. 



III. RESULTS 

In this section we will present the results. In the first part we will discuss the relaxation 
dynamics of the system at long times, i.e. the a— relaxation regime. In the second part we 
will focus on intermediate times, i.e. the /3— relaxation regime, and compare these results 
with the predictions of MCT. 

The quantity of main interest which will be studied in the present paper is the inter- 
mediate scattering function F{q,t) and its self part, Fs{q,t) These two space-time 



correlation functions can be defined most easily in terms of the density fiuctuations for the 
particles of type a 
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(^Pa (q, t ) = ^ exp (iq • (t) ) , (2) 

i=i 

where q is the wave-vector. The intermediate scattering function is then defined as 

F''^{q,t)= ] {6pM,t)Sp}{q,0)) . (3) 

I\Si + iVo 

Here (.) stands for the thermal average and we have assumed that the system is isotropic 
and hence F{q, t) depends only on the modulus q = |q|. From Eqs. and (0) it is clear that 
F{q,t) can be split into two parts, F°''^{q,t) = F^{q,t)6ai3 + Fl^^{q,t), where the functions 
F" and correspond to the diagonal and off diagonal elements in the double sum of F"'^: 

= ^E(e^P(^q- (r.W -r,(0)))) (4) 

and 

Ff{q,t) = EE-(exp(^q-(r,(t))-r,(0))) . (5) 

iVsi + iVo ;t = l 

Here the prime means that the term with k = j has to be left out. 

Note that Fs{q,t) and Fd{q,t) are the space- Fourier transforms of the self and distinct 



part of the van Hove correlation function |]2l|, thus quantities which are often studied in 



computer simulations of liquids. However, since in scattering experiments the latter func- 
tions are not directly accessible, in contrast to Fs{q,t) and Fii{q,t), and since also from a 
theoretical point of view the scattering functions are more interesting than the van Hove 
functions, it is appropriate to determine also them in simulations. We emphasize that to 
obtain the results of the present work we have not used the connection between the van 
Hove functions and the intermediate scattering functions, but calculated the latters directly 
from their definitions given in Eqs. (^) and (|^). For this we averaged over all q- vectors that 
had the same modulus and that were compatible with the size of the simulation box. 

In Fig. ID we show the time dependence of Fs{q,t) for the oxygen atoms for all tem- 
peratures investigated. The wave- vector is g = 1.7 A~^, the location of the first sharp 
diffraction peak in the structure factor, i.e. the length scale that corresponds to the typical 
distance between two adjacent tetrahedra in the network (We mention that for other 
wave- vectors the curves look qualitatively very similar.) From this graph we see that at 
high temperatures the curves show at very short times a crossover from a ballistic regime 
to a relaxation behavior which is basically exponential and that the correlation function 
decays to zero within one ps, i.e. very rapidly. To illustrate this exponential decay we have 
included in Fig. |l| an exponential function f{t) = 0.75 exp(— t/0.25ps), dashed line, which 
essentially coincides with the simulation curve at T = 6100 K for t > 0.1 ps. (We emphazise 
that due to the complex vibrational dynamics of the system at short times |[T^, one has 
to use an amplitude significantly smaller than 1.0) For intermediate and low temperatures 
one observes two additional regimes: Immediately after the ballistic regime, which last just 
up to ~ 0.02 ps, the motion of the ions shows a vibrational character, as can be inferred 
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from the presence of a dip in the correlator at around 0.2 ps. This type of dynamics is the 
dominant one for temperatures below the glass transition, but from the figure we recognize 
that it can be seen even at temperatures as high as 3580 K. From a microscopic point of 
view this dynamics corresponds to the rattling motion of the ions in the cage formed by 
their surrounding neighbors and the network. (Note that this cage is not rigid at all, since 
the particles that form the "walls" of this cage move themselves also.) Although every 
liquid will show some sort of vibrational dynamics, the one for Si02 is special in that it 
has also important contributions at quite low frequencies (1-2 THz) the origin of which has 
in recent years been the matter of a strong debate [rB,^|. For times which are somewhat 



longer than the ones for the vibrations, i.e. t > 0.2 ps, this type of motion is damped out 
and the particles slowly start to leave the mentioned cage. This dynamics, in the following 
called /5-process P,^, is thus the beginning of the a-process, i.e. the time regime in which 
the ions leave the cage completely and lead to the structural relaxation of the system. In 
the correlation functions it is seen as a plateau at intermediate times, the length of which 
rapidly increases with decreasing temperature thus pushing the a-relaxation to large times 
(see Fig. [^). MCT makes detailed predictions on the dynamics of the system on the time 
scale of the /3-relaxation and in the following we will check the validity of these in detail. 

From Fig. ^ one gets the impression that the shape of the curves in the a-relaxation 
regime does not depend on temperature, i.e. that the correlators obey the so-called time 
temperature superposition principle (TTSP). This means that a correlator 0(t,T) can be 
expressed by 

0(t,r) = 0(t/r(r)), (6) 

where t(T) is the typical time scale for the decay of the correlation function, i.e. the a- 
relaxation time, and is a master function. MCT predicts that there exists a so-called 
critical temperature Tc in the vicinity of which the TTSP is valid. Since this is one of the 
main predictions of the theory it is of interest to check whether or not it is valid for the 
present system. For this we have defined an a-relaxation time r by requiring that at time 
T the correlator has decayed to of its initial value. If the TTSP holds, a plot of the 
correlators at different temperatures versus t/riT) should give a master curve. In Fig. |^ 
we show this type of plot and we recognize that apparently the scaling does not work very 
satisfactory. To show this we have plotted in the inset of the figure an enlargement of the 
main figure at long rescaled times. From this inset we recognize that the fact that the TTSP 
is not fulfilled is not due to the noise in the data but rather a systematic trend in the curves. 
We have also made sure that this violation of the TTSP is not due to the way we have defined 



r, since a different definition gave qualitatively the same results Recall, however, that 
the TTSP is predicted to hold only slightly above Tc, which for the present system is around 
3330 K 1^. In fact, if the curves for the different temperatures are inspected carefully one 
finds that the ones in the temperature range 4700i^ > T > 34007^ (short dashed lines in 
Fig. 0) do indeed fall onto a master curve whereas the ones for high (T = 6100 K) and low 
(T < 3000 K) temperatures decay faster. Hence the curve at T = 2750 K (bold line in 
Fig. ^ is already very close to the exponentially decaying one at T = 6100 K (long dashed 
line in Fig. Hence we see that slightly above the correlation functions do indeed obey 
the TTSP, whereas it is violated if the temperature is too far below T^. This behavior is in 
qualitative agreement with the theory, since it is predicted that far above and far below Tc 
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the relaxation is of Debye type P,[T3 



We also mention that according to MCT the TTSP is supposed to hold also in the (3- 
relaxation regime, i.e. in the time window in which the correlation functions are close to the 
plateau. From the graph we see, however, that in this time regime the curves do not collapse 
at all. The reason for this is likely the fact that in this system the microscopic dynamics, 
i.e. the above mentioned low-frequency vibrations, are so pronounced that they completely 



obscure the TTSP in this time window |^|, in contrast to simple liquids whose vibrational 
dynamics does not have strong contributions at very low frequencies, i.e. around 1-2 THz. 
Below we will see, however, that certain predictions of MCT concerning the dynamics in the 
/^-relaxation regime are still valid, despite the presence of the strong vibrational dynamics. 

The results presented so far concerned the incoherent part of the intermediate scattering 
function of the oxygen atoms. For the case of the silicon atoms the time and temperature 
dependence is qualitatively similar, although for this species the TTSP at long times is 
fulfilled a bit better [|T^]. Since Fs{q,t) measures essentially the dynamics of a tagged 
particle within its cage and how it escapes this cage at long times, not much can be learned 
from this correlation function on the relative motion of the particles. Therefore it is of 
interest to study also the coherent intermediate scattering function Fd{q,t) which is defined 
in Eq. (^. In Fig. ^ we show the time dependence of the three partial F^^{q,t) for all 
temperatures investigated. The wave-vector is again 1.7 A~^, i.e. the location of the first 
sharp diffraction peak. Although from a qualitative point of view these correlation functions 
are similar to the incoherent ones, a closer inspection does show interesting differences. 
E.g. the coherent curves show the effects of the complex vibrational dynamics of the system 
much more pronounced in that, e.g., several maxima and minima can be seen between 0.02 ps 
and 1 ps, whereas the incoherent curves show only one local minimum (see Fig. |l[). We also 
mention that the modulus of the incoherent functions as well as and F^^ are always 
smaller than 1.0 since they are (normalized) autocorrelation functions. This is in contrast 
to Ff^'^ which, being not a autocorrelation function, can become larger than 1.0. A closer 
inspection of the curves for this correlator does indeed show that for times around 0.03 ps 
the curves are larger than 1.0, which is the reason why this correlator is very fiat at short 
times. 

From Fig. ^ one sees that in the a-relaxation regime the shape of the curves depends 
only weakly on temperature, i.e. that the TTSP is valid. That this is essentially the case is 
demonstrated in Fig. |, where we plot the coherent correlation functions for the case of 0-0 
versus the rescaled time t/r. We see that within the noise of the data the curves basically 
collapse onto one master curve. However, also in this case small systematic deviations from 
the TTSP, as discussed for the incoherent functions, can be found upon closer inspection 
of the curves. Thus this shows that also for these types of correlation functions the TTSP 
holds only for temperatures relatively close to T^. 

The intermediate scattering functions discussed so far can be directly measured in spin 
echo experiments or photon correlation experiments. However, in neutron or dynamic light 
scattering experiments only the time-Fourier transforms of these correlation functions are 
accessible. Therefore it is of interest to calculate these functions also, in order to see to 
what extent the various regimes seen in the time domain can be found in the frequency 
domain. Since the intermediate scattering functions have been measured over more than 
seven decades in time, the calculation of their Fourier transform is not a trivial matter. In 
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order to do this we have approximated each curve with a sphne under tension and calculated 



the Fourier transform of the spline by means of the Filon formula |2J]. Thus we obtained 
the dynamic scattering functions S{q,uj) and Ss{q,i^)- Since in the experiments one often 
measures the imaginary part of the susceptibility, x"{ly we have multiplied these functions 
with u to obtain cj) and respectively. 

In Fig. ^ we show Xs a function oi v = 2ituj for all temperatures investigated. We 
see that at 6100 K Xs shows only a single broad peak. The width of this peak at half 
maximum is approximately 1.5 decades which shows that it is not Debye-like, since this 
would correspond to a width of 1.14 decades. This increase is simply due to the fact that 
the broad peak at T = 6100 K is a sum of a microscopic and a relaxation peak, where the 
relaxation peak should be of Debye type in agreement with the exponential decay of the 
corresponding curve for Fs{q, t) in Fig. |1] for t > 0.1 ps. 

If the temperature is lowered the single peak splits up into two, i.e. into a microscopic 
peak at high frequencies and an a-peak at lower frequencies. With decreasing temperature 
the amplitude of the microscopic peak decreases but its location is independent of T. In 
contrast to this the a-peak quickly moves to lower frequencies and also its height increases 
slightly (to compensate the decrease in height of the microscopic peak). However, the shape 
of the peak seems to be essentially independent of T. At low temperatures we find a well 
defined plateau between the a-peak and the microscopic peak. A plot of log(Xs) versus 
log(z/) shows that at these low temperatures the whole spectrum can be described very well 
by a sum of a peak at microscopic frequencies and an a-peak with a high frequency wing that 
scales like z/~^, where the exponent y is ~ 0.6 and independent of T. Thus this result is in 
agreement with the von Schweidler law discussed below in the context of Fig. |1^ which gives 
an exponent 0.62. It is remarkable that this type of plot shows that for temperatures around 
3250 K, i.e. around T^, the susceptibility seems not to be just the sum of an a-peak and 
a microscopic peak, but that near the minimum between the two peaks x'li^) is enhanced. 
This means that there is an additional process in this frequency regime, the /3-process of 
MCT. We emphasize that this enhancement is only seen close to Tc, in agreement with the 
prediction of the theory. 

That the shape of the a-peak does indeed not depend on temperature, is demonstrated in 
Fig. Pd. Here we show (?) ^) I x'li.^^ J^max) as a function of zz/z/max, where z/max is the location 
of the maximum of the a-peak. From this plot we clearly see that within the accuracy of 
our data the curves for low temperatures fall on top of each other. Finally we mention that 
qualitatively the same results have been obtained for the case of the silicon atoms. 

Since the TTSP is valid to a high degree, we can determine the shape of the master 
curves as a function of the observable. We have done this for the case of the coherent and 
incoherent intermediate scattering function for several values of q. To determine this shape 
we proceeded as follows: For each correlator (^{t) we determined an a-relaxation time t' at 
2750 K by requiring that (f){t) = 0.1 [^. We then fitted (f){t) with a Kohlrausch- Williams- 
Watts (KWW) function 

0(t)=Aexp(-(t/rr) (7) 

where the amplitude A and the stretching parameter (3 are free fit parameters. This fit was 
done in the interval 0.02 < t/r' < 2 which for typical values of (3 corresponds to about 90% 
of the height of the a- relaxation (see Fig. We made sure that choosing a different interval 
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does not affect the results significantly and we estimate the systematic error bar of (3 due to 
this choice to be usually smaller than 0.03 |14| . In Fig. |^ we show the wave- vector dependence 
of j3 for the two incoherent functions as well as for two coherent ones. We see that for small 
and intermediate q the value of j3 is relatively large, in agreement with the observation that 



strong glass-formers do not show much stretching ||26|. However, for the large values of g, (3 
is significantly smaller than 1.0 since the curves for Fg decrease progressively with increasing 
q. This trend can be rationalized with the prediction of MCT that with increasing q 
the value of (3 should approach the value of the von Schweidler exponent fe, discussed in 
more detail below, which for our system is around 0.62, i.e. significantly smaller than 1.0. 
At g = 4 the value for (3{q) as determined from Fs for silicon and oxygen is around 
0.77, i.e. higher than the value predicted by the theory for large q. However, it has to be 
realized that g = 4 cannot be considered yet as "large" since at this wave- vector the 
static structure factor shows still significant oscillations 0. In a MCT calculation for a 
hard sphere system it has been found that only at very large wave-vectors the value of 13 
becomes close to the von Schweidler exponent, i.e. at q at which the static structure factor 
is essentially 1.0 [^]. In view of this result it is hence not surprising that the value of our 
/3 at the largest q is still significantly above b. 

The curves for the coherent functions oscillate around the ones for the incoherent ones. 
This oscillation is in phase with the static structure factor with peaks that correspond to the 
maxima in Saail)- Note that for large q the coherent functions are approximated well by 
the corresponding incoherent ones ||21[ and thus it is expected that the values of (3 become 
identical also. However, in the range of q for which we can reliably determine the value of 
13 this is not yet the case. This supports our argument that none of the [3{q) curves has yet 
converged to its asymptotic value for g — >■ oo. Note that qualitatively the same behavior 
for (3{q) has also recently been observed for water, where it was found that /5(g) converges 
to the von Schweidler exponent h for large g and that it shows oscillation in phase with the 
static structure for its coherent part PU| . 

We now address the temperature dependence of the a-relaxation time r. Since Si02 is a 
strong glass-former 0] one expects that r follows an Arrhenius law. In Ref. we have shown, 
however, that for the diffusion constant as well as the viscosity this law is found only at 
relatively low temperatures. For intermediate and high temperatures a significant deviation 
was found in that the temperature dependence was weaker than the one expected from an 
Arrhenius law. That this is the case for t(T) as well is demonstrated in Fig. ^ where we 
show the tt-relaxation time for Fs{q, t) for the oxygen atoms (the curves for the silicon atoms 
look qualitatively the same). The three curves that are shown correspond to g- values at the 
first sharp diffraction peak (g = 1.7 A~^), the location of the first minimum in the structure 
factor (g = 2.2 A~^), and the location of the main peak in S{q) (g = 2.8 A~^). We see that, 
apart from a g-dependent prefactor, the temperature dependence of the three relaxation 
times is very similar in that we find at low T an Arrhenius law and at intermediate and 
high temperatures a crossover to a weaker temperature dependence. The activation energies 
found at low T are between 5-5.5 eV, which compares well with the experimental values for 
the diffusion constant and the viscosity (4.7 eV and 5.33 eV, respectively [^,^). Thus we 
conclude that also with regard to this quantity the BKS model is quite reliable. 

In Ref. [0 we showed that the deviations from the Arrhenius law at intermediate and high 
temperatures are related to a change of the transport mechanism of the ions in that at low 
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T the atoms make independent hops whereas for high T the motion is much smoother and 
flow-hke. Such a change in mechanism is predicted by MCT around the critical temperature 
Tc which for our system has been estimated to be around 3330 K 0. The theory 

predicts that for temperatures a bit above the typical relaxation times increase like a 
power-law 

r(T)oc(T-T,)-^ (8) 

with an exponent 7 that is universal for the system, i.e. is independent of the observable 
(Si or O, value of g, etc.). Thus the bending of the curves in Fig. |^ can be rationalized by 
the crossover from the power-law behavior given by Eq. (||) to an Arrhenius behavior at low 
temperatures. 

MCT predicts that the exponent 7 is independent of the observable and that the diffusion 
constant D should be proportional to 1/r. In order to check the validity of this prediction 
we have calculated the product Dqt and plot it in the inset of Fig. ^ as a function of 1/T 
for the three relaxation times shown in the main figure. We see that although the product 
is basically constant at high temperatures, it increases continuously at intermediate and 
low temperatures. Since in this log-lin plot the curves at low T become a straight line we 
conclude that the activation energies for the diffusion constant and the relaxation times 
are not the same, a conclusion that is confirmed if one measures these energies directly 
(see below). Thus we conclude that the proportionality D oc does not hold around 
Tf. or lower temperatures. One reason for D 9^ might be the presence of dynamical 



heterogeneities I^Sj, i.e. that the cage that each particles sees changes significantly from 
particle to particle which has the effect that also the time to escape this cage depends 
strongly on the particle. It is then easy to see that in this case the (average) diffusion 
constant is not inversely proportional to the (average) relaxation time |Q. The observation 
that D is not strictly proportional to has in the past already been made in other systems 
for which the predictions of MCT hold extremely well not only on a qualitative but also on 
a quantitative level |3^. Thus, since the prediction of MCT that D oc r is only the result 



of a leading order calculation of the theory, the fact that the product Dt is not completely 
constant should not be taken as strong evidence that for the present system the theory is 
not applicable. 

Since for the present system the temperature dependence of the diffusion constant is not 
described so well by MCT, we investigate now the T-dependence of the relaxation times 
directly. In Ref. we determined the value of T^. to be around 3330 K. Therefore we plot in 
Fig. 1^ the a-relaxation time as a function of T — Tc, for those temperatures which are higher 
than Tc. From this plot we see that the data at high temperatures can indeed be described 
well with a power-law (bold straight lines), in agreement with the prediction of the theory 
(see Eq. (|D). The same conclusion holds for the diffusion constant of the oxygen atoms, 
as can be seen from the open circles in Fig. Similar results have also been obtained 
for the relaxation time of Fs{q,t) for silicon as well as the various coherent intermediate 



scattering functions [|14[. It should be noted, however, that this power law can be observed 
only in a temperature range in which r changes by about one decade. This is significantly 
less than the range that has been observed for simple liquids ||3^-p6|, water [^, molecular 



glassformers [ P7| , pS| , or polymeric systems Thus we conclude that for this system the 
corrections to the idealized version of the theory due to the hopping processes make the 
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observation of the power law rather difficult. 

As mentioned above, the value of the exponent 7 should be independent of the observable. 
In the inset of Fig. § we plot the value of 7 for the three different wave-vectors as well as 
the diffusion constant (point at g = 0). We see that whereas the value of 7 for intermediate 
values of q is around 2.4, the one for the diffusion constant is significantly smaller, 7 = 2.05. 
Thus again we find that this quantity behaves somewhat anomalous with respect to the 
ones at large g, and the reason for this is likely the one discussed above. Also included in 
the inset, horizontal dashed line, is the value of 7 as determined from a completely different 
approach, namely the analysis of the /3-relaxation regime (discussed below). We see that 
this estimate of 7 is in very good agreement with the values from the a-relaxation times, 
thus supporting the prediction of the theory that the exponent is not just a fit parameter 
but a quantity with a fundamental theoretical meaning. 

We now discuss the wave-vector dependence of the a-relaxation times in more detail. 
Since we have seen, see Fig. |^, that in the whole temperature range investigated the tempera- 
ture dependence of r is basically independent of g, it is sufficient to consider the g-dependence 
at only one temperature. This is done in Fig. |^ where we show the g-dependence of r' for the 
two incoherent as well as the three coherent functions. (Note, that for the reasons discussed 
in the context of Fig. ^ we show the g-dependence of r' instead of the one of r.) From 
the figure we see that the relaxation time for the incoherent functions, panel (a) and (c), 
decrease with increasing g. For a diffusive process, and hence also for small g, this decrease 
is given by 1/g^ [^. However, for the g-range considered here, i.e. g > 1.1 A"^, one would 
expect that the behavior of r is affected by the local structure. That this indeed the case 
is shown in Fig. |1^ where we plot r'g^ versus g. The product r'g^ can be interpreted as 
an inverse g-dependent diffusion constant. Thus, generally speaking, for g- values which 
correspond to length scales of nearest neighbor distances, r'g^ should be relatively large 
because then diffusion processes should be relatively slow due to the local order. As we see 
in Fig. [TO r'g^ exhibits an maximum around 2.7 A~^ for silicon as well as for oxygen which 
is close to the position of the main peak in the static structure factor corresponding to the 
length scale of a Si-0 bond. It is of course remarkable that no such feature is found at 
the position of the ffist sharp diffraction peak, at g = 1.7 A~^, although such a feature is 
probably buried under the broad peak around g = 2.7 (see Fig. |l^). We can rationalize 
this g-dependence from our previous observation [|^ that the elementary diffusion step for 
the oxygen atoms is due to the breaking of a Si-0 bond (and a slightly more complicated 
process for the silicon atoms) which shows that this length scale is of particular importance 
for the dynamics. We note that a similar behavior is found also in the MOT calculation for 
a hard sphere system PDI. In this case one obtains a peak in rg^ at a position which is close 



to the first peak in the static structure factor. 

In contrast to the monotonic decrease of the relaxation time for the incoherent functions, 
the one for the coherent functions shows an oscillatory behavior around the incoherent one, 
see Fig. These oscillations are in phase with the partial static structure factors, which 
have been included in the figures as well. Whenever the modulus of Sa/siq) is large, the 
corresponding r(g) is large also, a phenomenon that is very similar to the so-called de 
Gennes narrowing Note, however, that the de Gennes narrowing is a phenomenon 

observed at short times and thus cannot be used to explain the r(g) dependence on the 
time scale of the a-relaxation. Finally we mention that this type of dependence has also 
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been found in a MCT calculation for a hard sphere system and recently in a neutron 



scattering experiment of Toluene and m-Toluidine |41[ . 

Having discussed so far mainly the a-relaxation, we now turn our attention to the (3- 
relaxation, i.e. the time window in which the correlators are close to the plateau. MCT makes 
several interesting predictions about this relaxation regime and therefore we will check in 
the following to what extent they are valid for the present system. The first prediction of 
the theory is the so-called "factorization property", which says that close to the plateau 
every time correlation function (/>A(t) of an observable A is given by 

<PA{t) = <P\ + hAG{t), (9) 

where is the height of the plateau, Ha is a constant, and the whole time and temperature 
dependence is given by the function G'(t), which does not depend on the observable of A. If 
the factorization property holds one can see immediately that the ratio 

- Mt") - Mt') 

is independent of A if the times t, t', and t" are in the /3-relaxation regime, i.e. where Eq. 
holds. Note that this prediction of the theory does not only hold for the idealized version 
of the theory but also for the extended one, i.e. if the hopping processes mentioned in the 
Introduction are taken into account. Therefore it is reasonable to check the validity of 
Eq. (|^) not just above T^, but also below. This is done in Fig. |Tl] where we plot the ratio 



R{t) for a temperature above Tc, 4000K in panel (a), and below T^, 2750K in panel (b). The 
correlation functions shown are Fs{q,T) for silicon and oxygen at g = 1.7, 2.2, 2.8, 4.43, 
5.02 and 5.31 A~^. The times t" and t' are 0.4 ps and 1.6 ps for T = 4000 K and 11 ps and 
106 ps for T = 2750K, respectively. 

From Fig. pT| a we see that all the various curves collapse nicely onto a master curve at 
intermediate times, i.e. in the /3-relaxation regime. That this collapse is by no means trivial 
is seen at short and long times since there we find no collapse at all, thus showing that in 
general the shape of the correlation function does depend on the observable. Only in the 
/^-regime the correlators show a system universal time dependence, in agreement with the 
theory. 

The theory also predicts that with decreasing temperature the time window over which 
the /3-relaxation is observed should expand. That this is indeed the case is demonstrated in 
Fig. |TT]b, where we show the same correlation functions as in panel (a), but now at a lower 
temperature. From the figure we see that now the range over which the correlators fall onto 
a master function has indeed increased considerably, in agreement with the theory. Also 
included in this panel is an enlargement of R{t) at times where the correlator shows a dip, 
due to the boson peak mentioned above. From this inset we see that on this time scale the 
factorization property no longer holds since there the correlation functions are dominated 
by the vibrational dynamics for which Eq. is not valid. 

So far we have considered only the incoherent intermediate scattering functions. That the 



factorization property holds also for other time correlation functions is shown in Fig. |Tl|c 
where we plot, for T = 2750K, the time dependence of R(t) for Fs{q,t) for oxygen at 
q = 2.8 A~\ for Fd{q,t) for Si-Si at g = 1.7 and the functions aQ^(t) and aQ'*-*(t) which 



have been proposed by Roux et al. and which are defined as 



12 



a, 



Si(t) = / " ATTr'Gf{r,t)dr (11) 



and 



„SiO 




a^'"(t)= / 47rr2 G^'"(r,t)-1 G'^'"(r,0)-1 dr. (12) 



Here G^'(r, t) and (r, t) are the self and distinct part of the van Hove correlation func- 
tions 1^ which have been calculated and discussed in Ref. [0]. The values of and Tq'^ 
are 0.4 A and 2.35 A, respectively. Thus we see that the four time correlation functions 
considered are of very different nature and their overall relaxation behavior is certainly not 



the same. From Fig. p!l| c we see, however, that in the /3-relaxation regime the four functions 
fall nicely onto a master curve. Therefore we find that also for these correlators Eq. 
holds, which is hence strong evidence that for the present system the factorization property 
is a general property of a vary large class of time correlation functions. 

Having shown now that in the /5-relaxation regime the shape of a time correlation function 
4>Ait) is independent of the observable A considered, we continue now to study this shape 
in more detail. For this we make use of a further prediction of MCT, namely that the time 
dependence of (f)A{t) at long times is not arbitrary, but given by a sum of power laws of the 
form: 

Mt)=<pA-hAt' + hff' + ---. (13) 

Here 0^ is the height of the plateau, often also called the non-ergodicity parameter, and 
Ha (also called critical amplitude) and h^J^ are constants. According to MCT the value 
of the exponent b, also called von Schweidler exponent, is independent of A, i.e. it is a 
system-universal quantity. The first two terms on the right hand side are the so-called von 
Schweidler law. Note that Eq. (p!3[) is an expansion of the correlator in terms of t^, i.e. it is 
assumed that the third term is smaller than the second one. However, it has been suggested 
by the theory |^2[ that in order to determine reliably the coefficient Ha it is usually necessary 



to take into account these correction as well. In order to check whether the Ansatz (p!3D 
gives a good description of our correlation functions we have fitted them in the /^-relaxation 
regime with this functional form. The temperature we used was 2750 K, since this gives us 
the longest plateau and hence allows us to determine the fit parameters with the highest 



precision. The result of these fits are shown in Fig. |I2| for the case of the incoherent function 
for the oxygen atoms, panel (a), and the coherent function for the 0-0 correlation, panel 
(b). For each correlator (solid lines) two fits are shown: The one in which all three terms 
in Eq. (|13D were used (dotted lines) and the one in which the third term was neglected 
(dashed line). The difference between these two fits are seen at long times where the fits 
with the correction terms are able to describe the correlator for about one decade more than 
if these terms are not taken into account. That the quality of the fits improves significantly 
by considering the correction term oc t^** was first shown by Sciortino et al. for the case of 



water We emphazise that for these fits the value of h was kept fixed for all the different 
correlators, i.e. it was a global fit parameter. The value we obtained for h was 0.62 ±0.02. It 
can be seen from the figure that, if the correction terms are taken into account, the resulting 
fits are very good and describe the data over more than three decades in time. Therefore 
we conclude that this prediction of the theory, i.e. Eq. (|13D, holds for the present system. 
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A further interesting prediction of MCT is that there exists a one-to-one correspondence 
between the value of the von Schweidler exponent b from Eq. ([T3|) and the exponent 7 of 
the a-relaxation time in Eq. 0)0- According to the theory one can start with b and 
determine from 

r(l + 26) r(l-2a) 
the value of the parameter a and then use the relation 

7 = — + ^ (15) 
' 2a 2b ^ ^ 

to calculate 7. (Here r(a;) is the F-function.) If this is done for our value b = 0.62, one 
obtains 7 = 2.35. Since we have determined the value of 7, as defined by Eq. (|) also directly 
from the temperature dependence of the a-relaxation time, we can compare it now to this 
theoretical value. This is done in the inset of Fig. ^ and, as discussed already above, the 
theoretical value agrees well with the ones from t(T) for finite values of q. Thus we conclude 
that also this highly non-trivial prediction of the theory is satisfied quite accurately. 

In the last part of this section we discuss the wave- vector dependence of the parameters 
that where obtained from the fits of the correlators in the /3-regime to the functional form 
given by Eq. (|1^). For this we have made fits to the two incoherent and the three coherent 
intermediate scattering functions at T = 2750K for a wide range of wave-vectors. The de- 



pendence of the non-ergodicity parameters is shown in Fig. 0. For the case of the incoherent 
function, fg{q), we see that the g-dependence is Gaussian like. A fit with the functional form 
exp(— g^r^) gives a very good description of the data (see solid lines in panels (a) and (c)). 
For the localization length r^, i.e. the size of the cage, we obtain for silicon and oxygen the 
values 0.23 A and 0.29 A, respectively. This has to be compared with the typical nearest 
neighbor distance between silicon and oxygen, which is d = 1.6A [0,0. Thus we find that r<j 
is around 0.14d and 0.18d for silicon and oxygen, respectively. This size is somewhat larger 
than the one found in most solids which, in agreement with the Lindemann criterion P| 



have usually ^ O.ld — 0.14c/. The reason for this difference might lie in the fact that Si02 
is a very open network and thus the atoms can vibrate with a relatively large amplitude 
without approaching each other too much. 

The wave-vector dependence of the non-ergodicity parameter for the coherent functions 
shows oscillations which are in phase with the corresponding partial structure factor. For 
the case of Si-Si and 0-0 they oscillate around the corresponding and become basically 
indistinguishable from them at large q. These results are in qualitative agreement with the 
ones obtained for other systems [p9| , ^ , ^7H 39| which shows that from this point of view the 



present open network system does not behave differently as a simple liquid or a polymer. As 
already mentioned in the Introduction, it has very recently been shown that the g-dependence 
of the non-ergodicity parameters /'^(g) shown here can be calculated also reliably with the 
theory without any fit parameter. Since the agreement was also good at large wave- vectors, 
and since we see now that the coherent function oscillates around the incoherent one, we 
expect that the theory is indeed able to give a reliable estimate of the size of the cage also 
for this network forming system. 



The wave- vector dependence of the remaining parameters are shown in Fig. 14, In panel 



we show the critical amplitude hg and the first correction h^^^ for the two incoherent 



14 



functions. From this figure we recognize that from a quahtative point of view the curves are 
very similar to the ones calculated within the MCT for a system of hard spheres |2^. In par- 
ticular we find that h^^^ shows at small wave- vectors a negative minimum and subsequently a 
positive maximum at significantly larger q. In between it becomes zero at g = 1.98 A^^ and 
1.82 A^^ for the case of silicon and oxygen, respectively. This is close to the first peak in the 
structure factor, in agreement with the results for the hard sphere system p8|. This results 
shows that also from this point of view the theoretical calculation is qualitatively correct. 
Note that the fact that hf\q) = means that for this wave-vector the von Schweidler law, 
i.e. the second term in Eq. (|T^ can be seen over the broadest possible time-window, since 
the first correction to this law is zero. 

Also the g-dependence of h{q) and h^'^\q) for the collective quantities. Fig. |l^, is in 
qualitative agreement with the MCT results for a hard sphere system E.g. h^'^\q) 



shows at the location of the first peak a minimum and its value is negative. Thus we see 
that the theory is able to rationalize the behavior of this function, at least in a qualitative 
way. Of course it would be a more stringent test to do the MCT calculation directly for the 
silica system instead for hard spheres. However, this type of calculation is currently still 
very demanding and thus has to be left as a future project. 



IV. SUMMARY AND DISCUSSION 

In this paper we have presented the results of molecular dynamics computer simulations 
of a model for silica, the archetype of a strong glass former. The goal of this study was 
twofold: On the one hand to investigate in detail the relaxation dynamics of this system 
in the temperature regime where this dynamics is slow. On the other hand we wanted to 
check to what extent the mode-coupling theory of the glass transition is able to describe this 
dynamics at low temperatures. In these investigations we mainly focussed on the time-and 
temperature dependence of the various intermediate scattering functions for different wave- 
vectors. We found that these functions show nicely the slowing down of the system upon 
cooling in that they change their shape from a single exponential at high temperatures to 
a two-step relaxation at low temperatures. We emphasize that this two-step relaxation is 
already seen at temperatures which are about 30% above the melting temperature of the 
system, which shows that in order to show a slow dynamics the system does not have to be 
supercooled. 

In a previous paper it was demonstrated that this system shows a crossover of the 
dynamics from a flow-like motion at intermediate and high temperatures to a hopping-like 
motion at low temperatures |]^ . The temperature of this crossover can be identified with the 
critical temperature of MCT and is around 3330 K. We now find that the time-temperature 
superposition principle holds for temperatures around Tc and that time correlation functions 
are stretched. For significantly higher and lower temperatures, the TTSP is violated in that 
the correlators become more exponential-like. This observation is also in accord with the 
existence of the above mentioned crossover from a flow-like dynamics (i.e. collective motion 
which leads to a stretched exponential relaxation) to a dynamics in which hopping dominates 
(i.e. single particle dynamics which leads to an exponential relaxation). These results are in 
qualitative agreement with the prediction of the theory |13[. 

The temperature dependence of the a-relaxation times are qualitatively similar to the 
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one of the diffusion constant in that also they show a crossover from an Arrhenius law at low 
temperatures to a weaker temperature dependence at higher T. In this latter temperature 
range the T-dependence of r can be fitted well with a power law with the same critical 
temperature as found in Ref. [0]. However, the exponent 7 of this power law that we find 
for the intermediate scattering function is significantly higher than the one for the diffusion 
constant, in contradiction with the leading order prediction of the theory. Similar to the 
findings in simple liquids, the reason for this discrepancy might be the existence of dynamical 
heterogeneities, i.e. due to the fact that the temperature dependence of the diffusion constant 
is weaker than expected. 

According to the theory the time-temperature superposition principle should hold also in 
the /3-relaxation regime. We find that a plot of the correlation functions versus rescaled time 
does not lead to a collapse of the curves onto a master function in this time window. The 
reason for this is likely due to the presence of a very pronounced vibrational dynamics, related 
to the boson peak, that extends to relatively low frequency, i.e. influences the dynamics 
also at relatively long times. It has been shown before that the presence of a pronounced 
microscopic dynamics may disturb the /3-relaxation dynamics predicted by MCT so strongly 
that the asymptotic results cannot be observed anymore [^. Therefore it is not surprising 
that for the present system the TTSP does not hold in the /5-regime. (Note, that the reason 
that the TTSP is nevertheless observed in the a-regime is likely the fact that the latter 
regime is much less affected by the microscopic vibrations than the former since on the 
time scale of the a- relaxation the vibrations have finally been damped out.) However, by 
calculating the function R{t) from Eq. (|10D we can show that the time dependence of a wide 
class of correlation functions is the same for times in the /^-relaxation regime. In agreement 
with the theory this dependence is a sum of fractal power laws. We also find that the value 
of the von Schweidler exponent b in these laws obeys the connection proposed by MCT 
between b and the exponent 7 of the a-relaxation time, if one uses the 7 as determined from 
t{T) from the intermediate scattering function at finite wave- vectors. This result shows 
that neither 7 nor b are mere fit parameters but instead that they have a more fundamental 
theoretical meaning. 

By fitting the coherent and incoherent intermediate scattering functions in the j3- 
relaxation regime with the functional form proposed by MCT we have determined the 
wave-vector dependence of the various nonergodicity parameters, the critical amplitudes, 
and the correction terms. From a qualitative point of view these dependencies look very 
similar to the ones found in simple liquids or simple molecules [^,0,^,^. Since the lat- 
ter systems are all fragile glass formers we thus conclude that for these quantities there is 
no qualitative difference between fragile and strong glass formers. Furthermore it has to 
be mentioned that these dependencies are also in qualitative agreement with the ones pre- 
dicted by MCT for the simpler systems, thus giving evidence that the theory is indeed able 
to describe them correctly. 

In summary, we can say that MCT allows to understand many aspects of the dynamics 
of this strong glass former. Moreover, in Ref. [Q it was shown that the g-dependence of 
the Debye- Waller factor is predicted well by MCT. To what extent this is true also for the 
other observables studied in the present paper is not yet known and it is certainly important 
to check this in the future. It is of interest that we find that some of the predictions 
of the theory, such as the factorization property, hold also below Tc, i.e. in a temperature 
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regime where the hopping processes arc important. For the present system we find, however, 
that in the /9-relaxation regime the timc-tcmpcrature superposition principle docs not hold. 
Whether this is a real failure of the theory for strong glass formers or whether this is just 
a side effect of the strong vibrational dynamics at short times, can unfortunately not be 
decided. Hence it would be of great interest to investigate this question for a strong glass 
former which docs not show such a pronounced vibrational dynamics. 
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FIGURES 




FIG. 1. Time dependence of the incoherent intermediate scattering function for the oxygen 

atoms for all temperatures investigated. The wave- vector g is 1.7 A^^, the location of the first 
peak in the structure factor. The dashed line is a fit to the curve for 6100 K with a exponential 
function (see text for details). 
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FIG. 2. Incoherent intermediate scattering function for oxygen versus t/T{T), where r is the 
a-relaxation time at temperature T. Inset: Enlargement of the curves at long rescaled time. See 
text for the meaning of the different line styles. 
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FIG. 3. Time dependence of the coherent intermediate scattering function for all temperatures 
investigated. The wave-vector is 1.7 A~^, the location of the first peak in the structure factor, a) 
Si-Si, b) Si-0, c) 0-0. 
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FIG. 4. Coherent intermediate scattering function for the oxygen-oxygen correlation versus 
t/riT), where r is the a-relaxation time at temperature T. Inset: Enlargement of the curves at 
long rescaled time. See Fig. ^ for meaning of the dashed curves. 
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FIG. 5. Frequency dependence of the imaginary part of the dynamic susceptibility x'siQ^'^) 
for the oxygen atoms for all temperatures investigated, a) x'sil^^) versus frequency ly. b) 
Xsil' z^)/Xs(9, i^max) versus J^/fmax, where v^^x is the location of the a-peak in Xg- 




FIG. 6. Wave-vector dependence of the Kohlrausch- Williams- Watts stretching parameter for 
the silicon and oxygen incoherent intermediate scattering functions and the silicon-silicon and 
oxygen-oxygen coherent intermediate scattering functions. T = 2750 K. 
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FIG. 7. Temperature dependence of the a-relaxation time r' as determined from Fs{q,t) for 
oxygen for various values of q. The sohd straight hues are fits to the data at low temperature with 
an Arrhenius law with the activation energies given in the figure. Inset: Temperature dependence 
of the product Dqt', where Dq is the diffusion constant for the oxygen atoms. 
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FIG. 8. The a-relaxation time r', as determined from Fs{q,t) for oxygen, as a function of 
(T — Tc)/Tc, where Tc = 3330 K is the critical temperature from MCT as determined in Ref. 0. 
Also included is the inverse of the diffusion constant as determined in Ref. Q. The solid straight 
lines are fits with power laws of the form given by Eq. (|8|). Inset: Value of the exponent 7. The 
point at q = corresponds to the exponent for the diffusion constant. The horizontal dashed line 
is the value of 7 as determined from the dynamics in the /3-relaxation regime. 
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FIG. 9. Wave- vector dependence of the a-relaxation time r' at T = 2750 K as determined 
from the incoherent and coherent intermediate scattering functions (curves with symbols). Also 
included are the partial structure factors from Ref. |^ (curves with no symbols). 
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FIG. 10. The product T'{q)q^ as determined from the incoherent intermediate scattering 
functions as a function of g at T = 2750 K. 
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(a) 
T=4000K 
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FIG. 11. Time-dependence of the ratio R{t), as defined in Eq. (|lO|) , for various correlators, a) 
Fs{q,t) for silicon and oxygen at g = 1.7, 2.2, 2.8, 4.43, 5.02 and 5.31 A^^ T = 4000 K; b) the 
same set of correlators as in a) but now at T = 2750 K; c) four different correlators (see main text 
for description) at T = 2750 K. 
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FIG. 12. Time-dependence oi Fs{q,t) for oxygen, panel a), and F(i{q,t) for the oxygen-oxygen 
correlation, panel b), for various wave- vectors (solid lines). The dotted lines are fits the functional 
form given by Eq. ( [T3| ) . The dashed lines are fits with the functional form given by Eq. (|l^) without 
the last term. 
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FIG. 13. Wave- vector dependence of the nonergodicity parameters, a) fg{q) from Fs{q,t) 
for silicon and f^iq) from Fii{q,t) for the silicon-silicon correlation; b) f^^q) from F(i{q,t) for 
the silicon-oxygen correlation; c) fg{q) from Fs{q,t) for oxygen and f^i^q) from F(i{q,t) for the 
oxygen-oxygen correlation. The solid lines in a) and b) are a fit with a Gaussian to the incoherent 
data. 
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FIG. 14. Wave-vector dependence of the critical amplitude h{q) and the first correction term 
h^'^\q). a) Incoherent functions for silicon and oxygen; b) Coherent functions for the silicon-silicon 
and the oxygen-oxygen correlations. 
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